Gravitational Stability and Bulk 
Cosmology 

Nakia Carlevaro ^ and Giovanni Montani ^' ^' ^ 



'^Department of Physics, Polo Scientifico - Universitd degli Studi di Firenze, 
INFN - Section of Florence, Via G. Sansone, 1 (50019), Sesto Fiorentino (FI), Italy 

'^ICRA - International Center for Relativistic Astrophysics, 
c/o Dep. of Physics - "Sapienza" Universita di Roma 

Department of Physics - "Sapienza" Universitd di Roma, Piazza A. Moro, 5 (00185), Rome, Italy 

'^ENEA - C.R. Frascati (Department F.P.N.), Via Enrico Fermi, 45 (00044), Frascati (Rome), Italy 

ICRANet - C. C. Pescara, Piazzale della Repubblica, 10 (65100), Pescara, Italy 



nakia. carlevaroOicra. it montaniOicra. it 



Abstract: We present a discussion of the effects induced by bulk viscosity 
either on the very early Universe stabihty and on the dynamics associated to 
the extreme gravitational collapse of a gas cloud. In both cases the viscosity 
coefficient is related to the energy density p via a power-law of the form C = CoP^ 
(where Coj s = const.) and the behavior of the density contrast in analyzed. 

In the first case, matter filling the isotropic and homogeneous background is 
described by an ultra-relativistic equation of state. The analytic expression of the 
density contrast shows that its growth is suppressed forward in time as soon as Co 
overcomes a critical value. On the other hand, in such a regime, the asymptotic 
approach to the initial singularity admits an unstable collapsing picture. 

In the second case, we investigate the top-down fragmentation process of an 
uniform and spherically symmetric gas cloud within the framework of a Newtonian 
approach, including the negative pressure contribution associated to the bulk 
viscous phenomenology. In the extreme regime toward the singularity, we show 
that the density contrast associated to an adiabatic-like behavior of the gas (which 
is identified by a particular range of the politropic index) acquire, for sufficiently 
large viscous contributions, a vanishing behavior which prevents the formation of 
sub-structures. Such a feature is not present in the isothermal-like collapse. We 
also emphasize that in the adiabatic-like case bulk viscosity is also responsible for 
the appearance of a threshold scale (equivalent to a Jeans length) beyond which 
perturbations begin to increase. 

PACS: 98.80.-k, 95.30.Wi, 51.20.-hd 
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Introduction 



Both the extreme regime of a gravitational coUapse and the very early stages of the Universe 
evolution are characterized by a thermodynamics which can not be regarded as settled down 
into the equilibrium. At sufficiently high temperatures, cross sections of the micro-physical 
processes are no longer able to restore the thermodynamical equilibrium. Thus we meet 
stages where the expansion and collapse induce non-equilibrium phenomena. The average 
effect of having such kind of micro-physics results into dissipative processes appropriately 
described by the presence of bulk viscosity (. As shown in [TJ [H this kind of viscosity can 
be phenomenologically described by a function of the energy density p as C = CoP* where 
Co, s = const. 

In the first part of this work we investigate the effects that bulk viscosity has on the 
stability of the isotropic Universe [4], i.e. the dynamics of cosmological perturbations is 
analyzed when viscous phenomena affect the zeroth- and first-order evolution of the system. 
We consider a background corresponding to a Friedmann-Robertson- Walker model filled with 
ultra-relativistic viscous matter, whose coefhcient C corresponds to the choice 5 = 1/2 and 
then we develop a perturbation theory which generalizes the original analysis performed by 

Lifshitz [Hie]. 

As issue of our first analysis we find that two different dynamical regimes appear when 
viscosity is taken into account and the transition from one regime to the other one takes 
place when the parameter Co overcomes a given threshold value. However in both these 
stages of evolution the Universe results to be stable as it expands; the effect of increasing 
viscosity is that the density contrast begins to decrease with increasing time when is 
over the threshold. It follows that a real new feature arises with respect to the non-viscous 
(Co — 0) analysis when the collapsing point of view is addressed. In fact, if Co overcomes 
its critical value, the density contrast explodes asymptotically and the isotropic Universe 
results unstable approaching the initial singularity. Since a reliable estimation fixes the 
appearance of thermal bath into the equilibrium below temperatures O{10^^GeV) and this 
limit corresponds to the pre-inflationary age, our result supports the idea that an isotropic 
Universe outcomes only after a vacuum phase transition settled down. 

In the second part of this work, we study one of the most attractive challenge of relativistic 
cosmology concerning the research a self-consistent picture for processes of gravitational 
instability which connect the high isotropy of the cosmic microwaves background radiation 
with the striking dumpiness of the actual Universe. An interesting framework is provided by 
the top-down scheme of structure fragmentation which is based on the idea that perturbation 
scales, contained within a collapsing gas cloud, start to collapse (forming sub-structures) 
because their mass overcomes the decreasing Jeans value HI [3 [lU] of the background 
system. In a work by C. Hunter [TT], a specific model for a gas cloud fragmentation was 
addressed and the behavior of sub-scales density perturbations, outcoming in the extreme 
collapse, was analytically described. 

In this respect, we investigate how the above picture is modified by including, in the gas 
cloud dynamics, the presence of bulk viscosity effects [K]- We start describing the Lagrangian 
zeroth-order evolution by taking into account the force acting on the collapsing shell as a 
result of the negative pressure connected to the presence of bulk viscosity. We construct 
such an extension requiring that the asymptotic dynamics of the collapsing cloud is not 
qualitatively affected by the presence of viscosity: in this respect we can assume the viscous 
exponent as s = s/e. Then we face the Eulerian motion of the inhomogeneous perturbations 
living within the cloud. The resulting dynamics is treated in the asymptotic limit to the 
singularity. 

As a result, we show that the density contrast behaves, in the isothermal-like collapse {i.e. 
in correspondence of a politropic index 1 ^ 7 < ^/s), as in the non- viscous case Co = 0- On 
the other hand, the damping of perturbations increases monotonically as 7 runs from 4/3 to 



2 



5/3; in fact, for such adiabatic-like case, we see that density contrasts asymptotically vanish 
and no fragmentation processes take place within the cloud when the viscous corrections 
are sufficiently large. In particular, we observe the appearance of a threshold value for the 
scale of the collapsing perturbations depending on the values taken by the parameters Sind 
7 e (4/3, 5/3]; such a viscous effect corresponds to deal with an analogous of the Jeans length, 
above which perturbations are able to collapse. However such a threshold value does not 
ensure the diverging behavior of density contrasts which takes place, in turn, only when a 
second (greater) critical length is overcome. According to our analysis, if the viscous effects 
are sufficiently intense, the final system configuration is not a fragmented cloud as a cluster of 
sub-structures but simply a single object (a black hole, in the present case, because pressure 
forces are taken negligible). 



1 Bulk Viscosity Effect on the Early Universe 
1.1 General Remarks 

In order to describe the temporal evolution of the energy density small fiuctuation, we de- 
velop a perturbations theory on the Einstein equations. We limit our work to the study 
of space regions having small dimensions compared with the scale factor of the Universe 
a [13j . According to this approximation, we can consider a 3-dimensional Euclidean (time 
dependent) metric as spatial component of the background line element 

ds^ ^ dt^ - {dx^ + dy^ + dz^) . (1) 

In linear approximation, perturbed Einstein equations write as 



- 7:S;SR - 87rG<5r; , (2) 



where the term JT^ represents the perturbation of the energy-momentum tensor which de- 
scribes the properties of the matter involved in the cosmological collapse. The perturbations 
of the Ricci tensor SR'j^ can be written in terms of metric perturbations h'j^ = —Sg'^^, starting 
from the general expression for the perturbed curvature tensor. For convenience let us now 
introduce a new temporal variable 77, set by the relation dt = adrj, and use the symbol (') for 
its derivatives; we moreover impose, without loss of generality, that the synchronous reference 
system is still preserved in the perturbations scheme: /iqo = /iqq = 0- 

With the above assumptions, the perturbations of the Ricci tensor and the curvature scalar 
read: 

SR^o--i,k"-^h', (3a) 
SR^ = ^ih^''' -h-'p'^'), (3b) 



2a2 a3 2a3 ' ^^^^ 

SR =-^ihZ^'^-h-:')-\h"~%h' . (3d) 

By using these expressions, we are able to rewrite the left-hand side of Einstein eqs. ^ 
through the metric perturbations h'^. 
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Dynamical Representation of Perturbations Since we use an Euclidean background met- 
ric ([T]), we can expand the perturbations in plane waves of the form e"i "", where q is the 
dimensionless comoving wave vector being the physical one k = q/a- Here we investigate the 
gravitational stability properly described by the behavior of the energy density perturbations 
expressible only by a scalar function; in this sense we have to choose a scalar representation 
of the metric perturbations [6l[13]. Such a picture is made by the scalar harmonics Q = e"' '', 
from which the following tensor 

Qi^yiQ, = [|'5f-2M^]Q, (4) 

can be constructed. We can now express the time dependence of the gravitational perturba- 
tions through two functions X{i]), n{ri) and write the tensor in the form 

= A(7?)Pf + Kv)Qi , h = ^(77)Q . (5) 

1.2 Introduction to Viscous Cosmology 

The presence of dissipative processes within the Universe dynamics, as it is expected at 
temperatures above O{10^^GeV), can be expressed by an additional term in the standard 
ideal fluid energy-momentum tensor used in the original works by Lifshitz and Khalatnikov 
[SI [6^ . Using the conservation law T^. = we can express the viscous energy- momentum 
tensor as 

Tf^i, = {p + p)Uf,Ui, -pgf,^ , p = p-CuPp, (6) 

where p denotes the usual thermostatic pressure and is the bulk viscosity coefficient. In the 
homogeneous models this quantity depends only on time, and therefore we may consider it as 
a function of the Universe energy density p. According to literature developments [TJ [U [31 [T3] 
we assume that C depends on p via a power-law of the form 

C = CoP% (7) 

where Co ^ is a const, and s is a dimensionless parameter. Furthermore, in a comoving 
system the 4- velocity can be expressed as u'^ = ^/a, = and the viscous pressure p 
assumes the form 

p = p~3Cop'^'/a^ . (8) 

Let us now consider the earlier stages of a flat Universe corresponding to 77 ^ 1 and 
with the ultra-relativistic equation of state p — p/3. The Universe zeroth-order dynamics is 
described by the energy conservation equation and the Friedmann one, which are respectively 

p' + 3-{p + p)^0, ^ = ^^/3 7TGp. (9) 

In this analysis we assume s — 1/2 in order to deal with the maximum effect that bulk viscosity 
can have without dominating the Universe dynamics since it corresponds to a phenomeno- 
logical issue of perturbations to the thermodynamical equilibrium [131 [TC| . The solutions of 
the zeroth-order dynamics, for s = 1/2, assume the form 

p = Ca-(2+2") , a = ai 77^/" , = 1 - x Co , (10) 

being C an integration constant, x = ^/54nG and the parameter oi = {Scj^nCG/Sy^^^ . 
Since we consider an expanding Universe, the factor a must increase with positive power of 
the temporal variable {i.e. w > 0) thus we obtain the constraint ^ Co < Vx which ensures 
this feature. 
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1.3 Perturbation Theory in the Early Universe 

Let us now perturb the viscous energy-momentum tensor. Using the synchronous character 
of the perturbed metric we get the following expressions 

6T° ^ 5p , 5T^ = a{p + p) , (11) 
STi = <5f [-^'5p + C (<5uT, + h'/2a')] , (12) 

here = tj^ — 3Co sp*^^ a'/a^ and Vg is the sound speed of the fluid given by — Sp/5p. 

The presence of viscosity does not influence the expression of the Ricci tensor and its 
perturbations, thus we can still keep expressions ^ and use the perturbed form of the 
energy-momentum tensor to build up the equations which describe the dynamics of and 
6p. It is convenient to choose, as final equations, the ones obtained from the Einstein ones 
for a ^ j3 and for contraction over a and (3, which read respectively 
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MP + p) 

Furthermore, taking the 00-components of gravitational eqs., we get the expression of the 
density perturbations 

Substituting in eqs. ((131) , (HH) the zeroth-order solutions (fTU)) and the scalar representation 
of the metric perturbations (O, we can get, respectively, two equations for A, p which read 

^" + i^^'-T (^ + A^) = 0. (16) 



-H^(A + ,)(l + 3I^^) + g|^i^=0. 



(17) 



It is worth nothing that, among the solutions, there are some which can be removed by 
a simple transformation of the reference system, compatible with its synchronous character. 
They do not represent any real physical change in the metric and we must exclude them 
from the dynamics. Such fictitious solutions, in the ultra-relativistic limit, assume the form 
X — p ~ const, and A -|- ^ ~ l/if 13J. 

1.4 Behavior of the Density Contrast 

Let us now study the gravitational collapse dynamics of the primordial Universe near the 
initial Big-Bang in the limit rj 1. As in Lifshitz work [S] , we analyze the case of perturba- 
tions scale sufficiently large to use the approximation r^g <C 1. In our scheme eqs. ([TBI) and 
(|17p admit asymptotic analytic solutions for the functions A and p; in the leading order A 
takes the form 

^-^+C2, (18) 
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where Ci, C2 are two integration constants. Substituting this expression in eq. (|17p we get, 
in the same order of approximation, the behavior of the function ^ as 

f'-^ + C^' (19) 

where we have excluded the non-physical solutions introduced above. The constant Ci is 
given by the expression Ci = ^/b{3 — i/t^), A and B being constants having the form 

3 ^ ' iVC/Sa-^ - 3Co/w ' a;+" 

Let us now write the final form of perturbations pointing out their temporal dependence 
in the viscous Universe. The metric perturbations (O become 

hi = ^P^ + ^Qi + cAQ^ + P^) - (20) 

and, by p3|) and (fTU]). the density contrast 5 — ^p/p reads 

,5 = FiCirj'-'/- + C2rf + C^ri''-''- + Cirf-''-], (21) 

where C3 2,A/q^LjB, and F = u?Qq^ j'^. 

We have now to impose the conditions expressing the smallness of perturbations at a given 
initial time ryo- The inequalities /i^ ^ 1 and 5 <C 1 yield only two fundamental constraints 
for the integration constants: 

Ci « 77^/"-' , C2 « 1 , (22) 

for any w- value within the interval (0, 1]. Furthermore we find an additional condition 
which involves the wave number q and the integration constant C; in particular a rough 
estimation for a; < 1/3 of the inequalities C\ <^ 'r]o^'^~^ ^^'^ ^3 Vo^^^'^ yields the condition 
q < (GCr/o)"^/^" which ensures the smallness of the cosmological perturbations. 

Using the hypothesis 77 ^ 1 we can get the asymptotic form of the corrections to the cos- 
mological background. The exponents of the variable 77 can be positive or negative according 
to the value of the viscous parameter ix'(^o)- This behavior produces two different regimes of 
the density contrast evolution: 



Case ^ Co < Vsx Here perturbations increase forward in time. This behavior corre- 
sponds qualitatively to the same picture of the non- viscous Universe (obtained setting Co = 0) 
in which the expansion can not, nevertheless, imply the gravitational instability: if we con- 
sider the magnitude order rj ^ l/q, the constraints on Ci, C2 imply that S remains small 
even in the higher order of approximation. This behavior yields the gravitational stability of 
the primordial Universe. 



Case i/sx < Co < Vx In this regime the density contrast is suppressed behaving like 
a negative power of rj. When the density contrast results to be increasing, the presence of 
viscosity induces a damping of the perturbation evolution in the direction of the expanding 
Universe, so the cosmological stability is fortified since the leading rj powers are smaller than 
the non- viscous ones obtained setting Co = Cl- 
in this case the density fiuctuation decreases forward in time but the most interesting result 
is the instability which the isotropic and homogeneous Universe acquires in the direction of 
the collapse toward the Big-Bang. For Co > V^x the density contrast diverges approaching 
the cosmological singularity, i.e. for 77 — > 0. In this regime, scalar perturbations destroy 
asymptotically the primordial Universe symmetry. The dynamical implication of this issue is 
that an isotropic and homogeneous stage of the Universe can not be generated, from generic 
initial conditions, as far as the viscosity becomes smaller than the critical value Co = Vsx- 
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2 Spherically Symmetric Gas Cloud Fragmentation 



In this second part, we present an hydrodynamical analysis of a spherically symmetric gas 
cloud collapse. This model was firstly proposed in a work by C. Hunter [TT], where he 
supposed that a perfect fluid cloud becomes unstable with respect to its own gravitation and 
begin to condense. The collapsing cloud is assumed to be the dynamical background on which 
studying, in a Newtonian regime, the evolution of density perturbations generated on this 
basic flow. Such an analysis is suitable for the investigation of the cosmological structures 
formation in the top-down scheme [171 118] since it deals with the sub-structures temporal 
evolution compared with the basic flow of the gravitational collapse. 

2.1 Unperturbed dynamics 

The unperturbed flow is supposed to be homogeneous, spherically symmetric and initially 
at rest. Furthermore the gravitational forces are assumed to be very much greater than the 
pressure ones, which are therefore neglected in the zeroth-order analysis. In such an approach 
the gas results to be unstable, since there are no forces which can contrast the collapse, and 
the condensation starts immediately. 

The basic flow is governed by the Lagrangian motion equation of a spherically symmetric 
gas distribution which collapses under the only gravitational action. As in the previous 
approach (see eqs. ([6|), (O), in order to include bulk viscosity effects in the dynamics, we 
introduce the bulk pressure 

p = p-Cop'u^^, (23) 

where = (1,0) is the shell comoving 4- velocity. In the Newtonian limit we consider, the 
metric can be assumed as a flat Minkowskian one expressed in the usual spherical coordinates 
and the metric determinant g becomes g = —r^sin^d. In this case, for the 4-divergence wf^ 
we immediately obtain: uf^ — 2r/r . 

Considering the basic flow density a.s p — M/ {^/iirr^) and the pressure force acting on the 
collapsing shell of the form Fp = p Airr^ , the Lagrangian motion equation for a viscous fluid 
becomes now 

_ _GM_ _ C dr 

where C = 8tt(o {SM/AttY . The origin O is taken at the center of the gas, r is the radial 
distance, G the gravitational constant and M the mass of the gas inside a sphere of radius r. 
In what follows, we shall suppose that the gas was at a distance a from O in correspondence to 
the initial instant to; this distance a identifies a fiuid particle and will be used as a Lagrangian 
independent variable so r = r(a, t) and the following parametrization can be introduced 

r = acos2/3. (25) 

Here f3 = (3{t) is a time dependent function such that /3(to) — and (3{0) — 7r/2 since we 
choose the origin of time to have t = when r — and to takes negative values. 

Eq. ([M)) must be integrated to obtain the evolution of the radial velocity v — 9r/gt and the 
density p of the unperturbed flow. Let us now require that the viscosity does not influence 
the flnal form of the velocity which in the non- viscous Hunter analysis {i.e. Co = 0) is 
proportional to [H]. Substituting v — ^1 -Jr into eq. ¥2M) we see that, in correspondence 
to the choice s = s/e, it is again a solution as soon as i? = C — \/C^ + 2GM , where B 
assumes only negative values. Although this dynamics is analytically integrable only for the 
particular value s = s/e, the obtained behavior ~ i/v^ remains asymptotically (as r ^ 0) 
valid if the condition s < s/e is satisfied. 



7 



Using such a solution we are able to build an explicit form of the quantity /3 defined by 
(|25p . After standard manipulations we obtain the relation 

cos/3^ = 3A(-i) , (26) 

where A is defined to be A = — i3/2a^/^. 

It is more convenient to use an Eulerian representation of the flow field. To this end, 
once solving the well known Poisson equation for the gravitational potential (j), we obtain 
the unperturbed solutions describing the background motion; all these quantities take the 
explicit forms 

v==[i;,0,0], I) = -2r/3tan/3 , (27a) 
p^po cos-« /3 , (27b) 
(t)^ -2npoG{a^ ~r^/S)cos-^ P , (27c) 

where (') denotes the derivate with respect to time and the non-radial components of velocity 
must vanish since we are considering a spherical symmetry. 



2.2 First-order perturbation theory 

The zeroth-order motion of a viscous basic flow which collapses under the action of its own 
gravitation was discussed above. We shall now suppose that small disturbances appear on this 
field; the perturbations evolution can be described by the Eulerian equations in presence of 
dissipative processes which come out from the thermodynamical irreversibility of the collapse 
process and are due to the microphysics of non-equilibrium [T9\ . 

Perturbations are investigated in the Newtonian limit starting from the well known conti- 
nuity, Euler-Navier-Stokes and Poisson equations. In such a picture, we are now interested 
to study the effects of the thermostatic pressure on the perturbations evolution. We shall 
therefore consider terms due to the pressure p in the motion equations, which read 

/9 + V-(pv) = 0, (28a) 
v-l- (v- V) • V = -V0- VpVp + C/pV(V • v) , (28b) 
V^(t) = inGp . (28c) 

The gas is furthermore assumed to be barotropic, i.e. the pressure depends only by the 
background density by 

p = Kp\ (29) 

where k, 7 are constants and f ^ 7 ^/s. By this expression we are able to distinguish a set 
of different cases related to different values of the politropic index 7. The asymptotic value 
7=1 represents an isothermal behavior of the gas cloud; the case 7 = 5/3 describes, instead, 
an adiabatic behavior and it will be valid when changes are taking place so fast that no heat 
is transferred between elements of the gas. We can suppose that intermediate values of 7 will 
describe intermediate types of behavior between the isothermal and adiabatic ones. 

In this model, zeroth-order solutions (|27p are already verified since the pressure gradient, in 
the homogeneity hypothesis, vanishes and the pressure affects only the perturbative dynamics. 

Let us now investigate first-order fluctuations around the unperturbed solutions, i.e. we 
replace the perturbed quantities: (v -I- 6v), {p + 6p), {(j) + 6(j)) and {p + 6p) in eqs. ([25]) . 
Following the line of the Hunter work, equations for the perturbed quantities (5v, Sp and 
d(j) can be built. Skipping many technicality, we are able to manipulate the system of the 
motion eqs. and write down an unique equation which which governs the evolution of density 
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perturbations 6p in the viscous regime. Taking into account the parameterization ([25]) we 
get 



(cos^° pS'p - 6 sin/3 cos*^ PPSp) - AirGpo cos*^ (36 p 
= (vs^ cos^/3-6C/po sin/3cos"/3/3) D^dp+ 
+ C/po cos^^ ^ ^ 



(30) 



Here time differentiation is taken at some fixed comoving radial coordinate, Vg is the sound 
speed v1 = dp/ dp and is the Laplace operator as written in comoving spherical coordi- 
nates. 

In order to study the temporal evolution of density perturbations, we assume to expand 
5p in plane waves of the form 



5p{r, t) 5p (t) e 



— i k-] 



(31) 



where 1/k (with k — |k|) represents the initial length scale of the considered fluctuation. 
According to this assumption, we are able to write the asymptotic form of eq. ([30|) near the 
end of the collapse as (— t) — > 0. In this case, we can make the approximation sm/3 « 1 and 
the quantity cos/3 is given by (pS)) . 

The background motion equations were derived for a particular value of the viscosity pa- 
rameter s = s/e: substituting the basic flow density given by (|27b|l in the standard expression 
of the bulk viscosity ([7|), we obtain ( = CoPo''^ cos~^(3 . With these assumptions, eq. ([50]) 
now reads 
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y " 



i-t) 5q+ 



'16 A " 
9^2 



(3A)27-2/3 



2A 
3l 



where Vq = k7Pq ^ and the viscous parameter A is given by 

2.3 Asymptotic Analysis of The Density Contrast 



5q = 0, (32) 



(33) 



A complete solution of eq. ([5^ involves Bessel functions of flrst and second species J and Y 
respectively and it explicitly reads 

6g= CiGi{t) + C2G2{t) , (34) 

where Ci, C2 are integration constants and the functions Gi and G2 are deflned to be 

Gi(0 = (-i)-^+s^ J„ [q{-tY''~^] , (35) 
G2{t) = {-t)-'i+^ r4g(-t)V--7] ^ (36) 

having set the Bessel parameters n and q as 

n = [A^ - 2AA + A^ + levrGpo]* / (6^(4/3 - 7)) , 



We now proceed, in order to study the asymptotic evolution of the solution , analyzing 
the cases 1 ^ 7 < Vs and ^/s < 7 ^ separately, since Bessel functions have diff^erent limits 
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connected to the magnitude of their argument. In the asymptotic hmit to the singularity, the 
isothermal-hke case is characterized by a positive time exponent inside Bessel functions so 
g(— t)*/^"''' <C 1, on the other hand, in the adiabatic-hke behavior, we obtain q{—ty^^~'^ ^ 1. 
For an argument mush less than unity Bessel functions J and Y behave like a power-law of the 
form Jn{x) ~ x+", Y„(x) ~ x~", for a; <C 1. And in the adiabatic-like case argument becomes 
much gather than unity and they assume an oscillating behavior like Jn{x) ~ x~^/^cos(a;), 

Y„{x) ~ x~^^'^sin{x), for X ^ 1. 



Isothermal-like Case In this first regime, an asymptotic form of functions G can be found 
as follow 

Gif = ci (-i)-f +B^±(|-^)" (38) 

where ci and C2 are constants quantities. The condition which implies the density pertur- 
bations collapse is that at least one of G functions diverges as {—t) 0. An analysis of 
time exponents yields that Gi diverges if A < 7A — 2nGpo/'iA but, on the other hand, G2 
is always divergent for all A. These results imply that, in the isothermal case, perturbations 
always condense. 

Let us now compare this collapse with the basic flow one; the background density evolves 
like cos'^P (see JSEl) that is, using ^ 

p-{~t)-^. (39) 

In the non-dissipative case (A — 0) perturbations grow more rapidly with respect to the 
background density involving the fragmentation of the basic flow independently on the value 
of 7; in presence of viscosity the density contrast assumes the asymptotic form 

SiSO _ j._^^-i + 5^-^y[A2-2AA+A^ + 167rGpo] _ 

Here the exponent is always negative and it does not depend on 7, this implies that (5^'^'^ 
diverges as the singularity is approached and real sub-structures are formed involving the 
basic flow fragmentation. This issue means that the viscous forces do not have enough 
strength to contrast an isothermal perturbations collapse in order to form of an unique 
structure. 



Adiabatic-like Case For 4/3 < 7 ^ 5/3^ j and Y assume an oscillating behavior. In this 
regime functions G read 

G,^r = lZH~tf"-''] (-i)*-^+^ , (41) 

where Ci_2 are constants. Following the isothermal approach, we shall now analyze the time 
power-law exponent in order to determine the collapse conditions. G functions diverge, 
involving perturbations condensation, if the parameter A is less than a threshold value: this 
condition reads A < 17A — 3A7 (for a given value of the index 7). Expressing A in function of 
the wave number (j33p . wc outline, for a fixed viscous parameter a constraint on k which 
is similar to the condition appearing in the Jeans model [7], [S]. The threshold value for the 
wave number is given by the relation 

kl = {nA-2,^A)pf IC^ (42) 

and therefore the condition for the density perturbations collapse, i.e. 5p^^^ 00, reads 
k < kc, recalling that, in the Jeans model for a static background, the condition for the 
collapse is fc < A:j = [inGpo / v'iY^^- It is to be remarked that, in absence of viscosity 
(Co — 0), expression (j42|) diverges implying that all perturbations scales can be conducted 
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to the collapse. On the other hand, if we consider perturbations of fixed wave number, they 
asymptotically decrease as (— t) ^ for A > 17 A — 3Aj. Thus for each k there is a value 
of the bulk viscosity coefficient over which the dissipative forces contrast the formation of 
sub-structures. 

If fc < kc, perturbations oscillate with ever increasing frequency and amplitude. For a 
non-zero viscosity coefficient, the density contrast evolves like 

A study of the time exponent yields a new threshold value. If A < 5A~3A'j, i.e. the viscosity 
is enough small, sub-structures form; on the other hand, when the parameter Co i or the wave 
number k, provides a A-term overcoming this value, the perturbations collapse is so much 
contrasted that no fragmentation process occurs. In other words, if A > 5^ — SA'f we get 
gADB _^ j_g^ £qj, ^ given 7 there is a viscous coefficient enough large ables to prevents 
the sub-structures formation. 

It is remarkable that in the pure adiabatic case, 7 = ^/s, dissipative processes, of any 
magnitude order, contrast the fragmentation because, while the Jeans-like length survives, 
the threshold value for sub-structures formation approaches infinity. We can conclude that, 
in the case Vs < 7 ^ ^/s^ the fragmentation in the top-down scheme is deeply unfavored by 
the presence of bulk viscosity which strongly contrasts the density perturbations collapse. 

3 Concluding remarks 

To complete the analysis of the gravitational stability in presence of dissipative affects, we 
now point out some relevant feature about the validity of our picture concerning both the 
relativistic analysis and the Newtonian fragmentation: 

Shear Viscosity In our approach, the hypothesis of the Universe isotropy and the fact 
that the shell results comoving with the collapsing background, imply that there are no 
displacements between parts of fluid with respect to ones other in both cases. Dissipative 
processes are therefore related only to thermodynamical properties of the fluid compression 
and can be phenomenologically described by the presence of bulk viscosity. Furthermore, 
in this model we are able to neglect the so-called shear viscosity since it is connected with 
processes of relative motion among different parts of the fluid. 

Early Universe Dynamics In the fist Universe analysis, we fix the value s = 1/2 in or- 
der to deal with the maximum effect that bulk viscosity can have without dominating the 
dynamics. In fact, the notion of this kind of viscosity corresponds to a phenomenological 
issue of perturbations to the thermodynamical equilibrium. In this sense, we remark that, 
if s > 1/2, the dissipative effects become dominant and non-pcrturbative. Moreover, if we 
assume the viscous parameter s < 1/2, the dynamics of the early Universe is characterized by 
an expansion via a power-law a{t) ^ t^^'^'^ starting from a perfect fluid Friedmann singularity 
at i = (here 7 is identify by the relation p = (7 — l)e)- After this first stage of evolu- 
tion, where viscosity does not affect at all the dynamics, the Universe infiates in the limit 
t ^ 00 {i.e. out of our approximation scheme) to a viscous deSitter solution characterized 
by a{t) - e"'>\ Hq = = yaiCoV^/lV^^'^'^''' Since, in this work, we deal with 

the asymptotic limit t — > 0, we only treat the case s = 1/2 in order to quantitatively include 
dissipative effects in the primordial dynamics. 

Zeroth-order Cloud Dynamics We now clarify why the choice s = s/e is appropriate to 
a consistent treatment of the asymptotic viscous collapse. We start by observing that bulk 
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viscous effects can be treated in a predictive way only if they behave as small corrections to 
the thermodynamical system. To this respect we have to require that the asymptotic collapse 
is yet appropriately described by the non viscous background flow. Taking into account eq. 
(|24|) it is easy to infer that in the asymptotic limit as r ^ 0, the non viscous behavior 
V ^ is preserved only if s ^ s/e; in fact, in correspondence to this restriction, the viscous 
correction, behaving like 0{r~^^^^/^), is negligible with respect to the leading order 0{r~^) 
when the singularity is approached. On the other hand, it is immediate to verify that, as 
(— i) —^ 0, the viscous terms in Sp and in Sg respectively are negligible and the perturbation 
dynamics matches asymptotically the non-viscous Hunter result if s < ^/e. Matching 
together the above considerations for the zeroth- and first-order respectively, we see that 
s — 5/6 is the only physical value which docs not affect the background dynamics but makes 
important the viscous corrections in the asymptotic behavior of the density contrasts. 

Validity of the Newtonian Approximation Since the analysis of the cloud collapse ad- 
dresses Newtonian dynamics while the cloud approaches the extreme collapse, it is relevant 
to precise the conditions which ensure the validity of such a scheme. The request that the 
shell corresponding to the radial coordinate r lives in the Newtonian paradigm leads to im- 
pose that it remains greater than its own Schwarzschild Radius, r{t) ^ 2GM{a) , where 
M{a) = po (|7ra^). About the dynamics of a physical perturbations scale I — (27r/fc)cos/3^ 
(here cos0^ plays the same role of a cosmic scale factor), its Newtonian evolution is ensured 
by the linear behavior, as soon as, the Schwarzschild condition for the background holds. 
More precisely a perturbations scale is Newtonian if its size is much smaller than the typical 
space-time curvature length, but for a weak gravitational field this requirement must have 
no physical relevance. To explicit such a condition, we require that the physical perturba- 
tions scale is much greater than its own Schwarzschild Radius, which leads to the inequality 

1 /2 

k ^ x(— where x — [f (27i')'^Gpo(3v4)^^/'^] . Such a condition tells us which modes 
are Newtonian within the shell. 
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